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Abstract 

We consider the effect of regularization by interval truncation on the spectrum of a singular 
non-selfadjoint Sturm-Liouville operator. We present results on spectral inclusion and spectral 
exactness for the cases where the singularity is in Sims Case II or Sims Case III. For Sims Case 
I we present a test for spectral inexactness, which can be used to detect when the interval trun- 
cation process is generating spurious eigenvalues. Numerical results illustrate the effectiveness 
of this test. 
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1 Introduction 

Over the last 30 years there has been considerable interest in numerical solution of singular Sturm- 
Liouville problems, and in particular in the development of automatic software for such problems: 
see, e.g., Bailey, Gordon and Shampine M, Bailey, Garbow, Kaper and Zettl j|, Fulton and Pruess 
[|| and Marietta and Pryce Q. The software described in these papers usually uses an interval 
truncation procedure to regularize problems posed either on infinite intervals, or on finite intervals 
with singular behaviour of the coefficients near at least one of the endpoints. Rigorous mathematical 
justification of the validity of the interval truncation process, however, did not appear (except 
for special cases) until the paper of Bailey, Everitt, Weidmann and Zettl H in 1993, which uses 
fundamental ideas from Reed and Simon ]l5| to develop conditions under which the spectra of 
a sequence of regularized problems can (a) provide approximations to the whole spectrum of the 
original singular eigenvalue problem (spectral inclusion); (b) not yield approximations to any points 
which are not in the spectrum of the original singular problem (spectral exactness). All of this work 
is for self adjoint problems only. 

Non-selfadjoint singular problems are also very important. They arise when the complex scaling 
method is used to find resonances of a selfadjoint problem (for a review see [jl2)) and also, more 
classically, in the study of hydrodynamic stability, where the spectra of the Orr-Sommerfeld and 
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related equations are often studied over infinite intervals. Recent applications to the study of large 
disturbances in water waves are described by Chamberlain and Porter |1. 

It is well known that the spectra of non-selfadjoint operators can be pathologically sensitive to 
perturbation of the operator. Matrix examples of such sensitivity are provided in the classic text 
of Wilkinson ^0|. For a recent study in the context of non-selfadjoint Sturm-Liouville operators 
see Davies ||, and for a study in the context of general operators via pseudospectra see Trefethen 
jl7|. Given this sensitivity, it seems important to ask: under what conditions can one expect the 
regularization process used for selfadjoint singular Sturm-Liouville operators to be successful for non- 
sclfadjoint Sturm-Liouville operators? In particular, can one recover results on spectral inclusion and 
spectral exactness? If not, might one at least be able to recover results on pseudospectral inclusion 
and pseudospectral exactness, or develop a-posteriori tests for spectral exactness? 

We seek to answer these questions in this paper, for singular second order non-selfadjoint Sturm- 
Liouville problems. 

For selfadjoint problems a singular endpoint is either of limit point or of limit circle type. This 
is the Titchmarsh-Weyl theory and may be developed either using methods of complex analysis (see 
Titchmarsh Jl8|) or using the theory of deficiency indices for symmetric operators on Hilbert spaces 
(see, e.g., Dunford and Schwartz H). The analogous theory for non-selfadjoint problems is due to 
Sims and to Brown, Evans, McCormack and Plum [[|. It is based on the Titchmarsh approach 
to the selfadjoint case, and will be very important in this paper. The other ingredient which we 
shall find useful is a non-selfadjoint analogue of the results of Reed and Simon on spectral inclusion 
and spectral exactness |l5|, theorems VIII. 23- VIII. 25], for which we shall use results from Harrabi 
[0 and Kato g p. 208]. 

2 A review of the Sims Classification 

The problem which we consider concerns the spectral behavior of 

M[y] = -[-(py)' +qy] on [a, 6), (1) 

where as usual 

(i) w > 0, p ^ a.e. on [a, b) and w,l/p G Lj oc [a,b); 

(ii) p, q are complex- valued, q G Lj oc [a,b) and 

Q = coi-^r +rp(x) :xe[a,b), < r < oo\ ^ C. (2) 
[w{x) J 

These assumptions imply that a is a regular point of ([!]) and we shall assume that b is a singular 
point. By this we mean that either b = +oo or that J a (w+ I Q \)dx — oo. Since we are assuming 
that Q does not occupy all of C, it is known that its complement has either one or two connected 
components. For Ao G C\Q we denote by K = K(Xo) the nearest point in Q to Ao and by L the 
tangent to Q at K and arrange by translation and rotation through an angle rj for L to coincide with 
the imaginary axis while Ao and Q are contained in the new left and right half planes respectively. 
That is, for all x 6 [a, b) and r £ (0, oo), we require by choice of K and t] that 

$t[{rp(x) + 44 - K}e ir >] > (3) 
w{x) 

and 

5R[(A - K)e 1 ^) < 0. (4) 
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The set of all such admissible pairs (rj, K) we call S and we also define 

A V ,K = {A G C : R[(A - K)e ir >\ < 0}. 

In order to obtain from (|l|) a well posed eigenvalue problem we need to introduce boundary conditions 
at a and possibly at b. The conditions at a will be given in the form 

y(a) cos a + py'(a) sin a = 0, (5) 

where the parameter a, which may be complex, will be subject to the condition 



Mucosa sin a] < 0. (6) 

This gives rise to a set S(a) which is defined as the subset of S in which (ph holds. We note that 
a = and a = ir/2 correspond to Dirichlet and Neumann boundary conditions resp ectively. 

When p — 1 and q is real the classical theory of Weyl jl9| and Titchmarsh Jll| shows that if 9 
and <f> are linearly independent solutions of (|l|) which satisfy 

4>{a, A) = sin a, 9(a, A) = cos a, 
p<f> (a, A) = — cos a, p9 (a, A) = sin a, (7) 

where a is now real, then there is a complex number m(A), a function of the strictly complex variable 
A, such that 

xjj = 9 + m<p (8) 

lies in L 2 [a, 6). When, up to constant multiples, ip is the only solution of the differential equation 
which lies in L^[a, 6), we say that (Q) is in the limit point case at b. If however both 9 and (p lie 
in L^Jfl, b) then we say that (|l|) is in the limit circle case at b. In this case an additional boundary 
condition at b is needed in order to make into a well-posed eigenvalue problem. There is a 

one to one correspondence between this additional boundary condition and the choice of function 
m(-) in (|§|), in the sense that with an appropriate choice of boundary condition there exists a unique 
function m(-) such that eqn. (|J) defines a solution ip of the differential equation satisfying the 
boundary condition at x = b, while with an allowed choice of m(-) the function ip defined in (|J) can 
itself be used, for appropriate A, to define the boundary condition at x — b in the form [y, ip](b) = 0, 
where [/, g] := p{fg' — f'g) denotes the Wronsian of two functions / and g. It is known that the 
classification of limit point or limit circle is independent of the strictly complex parameter A. The 
terminology of limit point or limit circle owes its origin to the method used to establish the existence 
and possible uniqueness of ip in (||) . It may be shown that the spectral points of any realisation of 
(|j) as an operator in l? w [a, b) may be charecterised by the behaviour in the limit as 3A — * of the 
function m(A) associated with the boundary conditions defining the domain of the realisation. 

Many of these notions may be carried over to the case when p, q and a are complex. In a 



seminal paper Sims |16 shows that when p = w = 1 and 3g € C_, where C_ denotes the strictly 
lower complex plane, then the limit point / limit circle classification of Weyl now gets replaced by a 
threefold classification. We shall discuss this in the more general setting of Q which only requires 
@ an d (f§) to hold. Using a nesting circle method based on that of both Weyl and Sims, Brown 
et al. prove the following theorem. 

Theorem 2.1 /Q/ For A G A^ ^, {rj,K) G S(a) the following distinct cases are possible, the first 
two being sub-cases of the limit point case: 

• Case I : there exists a unique solution of ^) satisfying 

W V {P I V I 2 +(V - Kw) | y \ 2 }}dx + | y | 2 wdx < oo; (9) 

J a 

and this is the only solution satisfying y G L^[a,b); 
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• Case II : there exists a unique solution of ^) satisfying ^), but all solutions of (Qj lie in 
L 2 Ja,b); 

• Case III: all solutions of ^) lie in L^[a,b) and satisfy (^j. 

It may also be shown that the classification is independent of A in the sense that 

(i) if all solutions of (jl[) satisfy (||) for some A G ^-ri,K Case III) then all solutions 
°f (OP sa >ti s fy (H^ f or a M A G C; 

(ii) if all solutions of ^) lie in L 2 w \a, b) for some A G C then all solutions of (Q) satisfy 
y G L 2 w [a,b) for all A G C. 

It is interesting to examine the case when p is real and non-negative. In this case for some rj G [—§,§] 
and K € C let 

9 K>n (x) = K[e'" ) (g(a;) - Kw{x))} > a.e. x G (a, b). (10) 

Then the condition (||) in the Sims characterisation of (|l|) in Theorem [2 .l| for A G A n .x, (v-> K) e S(a), 
becomes 

cosr] p\y' | 2 dx + Kr) (x) | y(x) 2 dx + / | y(x) | 2 w(x)dx < oo. (11) 

J a J a J a 

In this case the remark on the independence of the classification can be extended to the following: 



(i) if for some A G C all the solutions of (|lj) satisfy (|ll|), then for all A G C all solutions of 
(0) satisfy ©; 



(ii) if for some A G C all the solutions of (Q) satisfy one of 

,b 

cos r] / p I y' | 2 dx < oo, (12) 

J a 

,6 

/ 6 Kv y | 2 dx < oo, (13) 

J a 

then the same applies for all A G C. 

We remark that Sim's analysis is the special case of the above when rj ~ 7r/2, K = 0. This restriction 
overlooks the interesting features present in ( p"l| ) when r/ G (— 5, ^ ), namely, that the classification 



in Theorem 2.1 involves a weighted Sobolev space as well as L^[a,b). The paper JlJ also examines 
the analytic behaviour of m(A) and the connection between this and the spectrum of M, an operator 
realisation of Ad in L^a,^). This is sumarised in the following theorems, in which m denotes the 
unique function such that (^) defines a solution of ([!]) which either (i) lies in -^[a, b) (Sims Case I) 
or (ii) lies in L^[a,b) and satisfies the additional boundary condition at x = b (Sims Cases II and 
III). 

Theorem 2.2 (Theorem 4-7, ^) In Cases II and III, A is a pole of m of order s if and only if A 
is an eigenvalue of M of algebraic multiplicity s. 

Theorem 2.3 (Theorem 4.13, Suppose that (||) is in Case I. Define 

Q(a) = p| (C\A n>K ), 
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and let Q c (a) denote the set Q(a) when the underlying interval is [c,b) rather than [a,b). Define 

Q c := co { +rp{x) : x G [c, b), r G (0,oo)l , 

{w(x) J 

Qb ■= ^c.e(a,b)Qc, Qb{ct) = ^c£(a,b)Qc(a), 

Then m(X) is defined throughout C\Q(a) and has a meromorphic extension to C\Qb(ot), with poles 
only in Q(a)\Qb(a). In addition X is a pole of m(X) if and only if X is an eigenvalue of M for 
X^Q b {a). 



3 Tests for spectral inclusion and spectral exactness 



In this section we prove a simple theorem (Theorem |3.1[ ) which allows us to test a convergent 
sequence of eigenvalue approximations obtained from a sequence of truncated interval problems, in 
order to determine whether or not the limit of the sequence is truly an eigenvalue of our original 



problem. We also prove two additional results (Theorem |3.3| and Theorem 3.4) which give methods 
for determining whether or not the hypotheses of Theorem |3.l| are satisfied for a given problem. 
The second of these, Theorem [O], extends a convergence result in || from the complement of the 
numerical range of our singular operator into a set which is typically much larger. 



Finally, we show how Theorem 3.3 and Theorem 3.4 allow us to develop a test for spectral inclu- 
sion, to ensure that there will be no eigenvalues of the original problem which remain unapproximated 
by the truncation process. 

3.1 Spectral Exactness 

We denote by m(-) the m-function developed in section @; in Sims cases II and III, this corresponds 



to a particular choice of boundary condition at x = b. From Theorems 2.2 and 2.3 we know that 
the poles of m are the eigenvalues of a realization M of the differential operator A4 subject to a 
boundary condition of the form 

(cosa)y(a) + (sina)py'(a) = (14) 

(and possibly an additional boundary condition at x = b). 

We denote by L the operator, and by £(■) the Titchmarsh-Weyl function, when the boundary 
condition at x = a is changed to 

(sin a) y (a) — (cosa)py'(a) = 0, (15) 

but the boundary conditions at x = b (where applicable) are left unchanged, the same as those for 
M. 

The functions t and m are related by the identity 

m(X)l(X) = -1 

(see [|§ eqn. 5.17]). 

Let (6„)„ e N be a sequence such that b n f b as n /* oo. Following j|, §2] we may construct 
a sequence M n of regular operators defined on the intervals [a, b n ]. These operators M n are still 
defined by M n y = M.y on their domains: it is the boundary conditions defining these domains which 
are of interest. At x — a we keep the boundary condition (|l4|). At x = b n we impose a boundary 
condition of the form 

y(b n ) cos I3 n + py' (b n ) sin f3 n = 0. (16) 
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The Titchmarsh-Weyl function m n associated with M n is then given in terms of the solutions 6 and 
4> of (0) by 

m " 1 ' cf>(b n ,\) cot (3 n + P (f>'(b n ,\y 

[g eq. (2.6)]. We shall examine in Lemma |3.2| , Theorem |3.3| and Theorem [3.4| below conditions on 
the P n which ensure that, for A in certain regions of C, 



m(A) = lim m n (A). 



(17) 



Let L n be a sequence of regular operators defined on the intervals [a,6„] with boundary condition 
( [l5|) at x — a and with the same boundary conditions (]l^) as the M n at x = b n , so that the 
associated Titchmarsh-Weyl functions l n {') satisfy 

m n (X)£ n (X) = -1. 

By analogy with (|l7|) we shall assume that in some appropriate regions of C, 

£(A) = lim £„(A). (18) 



Theorem 3.1 (Test for spectral inexactness) Suppose that /i 6 C has the following properties: 

1. fi does not lie in the spectrum of M; 

2. m n —> m and £„—>■£ uniformly on any compact annulus of sufficiently small outer radius 
surrounding fi. 

Then there are only two possibilities: 

(a) there exists a neighbourhood N of n and N € N such that no eigenvalue of M n lies in AT for 

any n > N; 

(b) there exists a monotone increasing sequence (rij)j^ of positive integers, and two associated 

sequences (Xj)j^, (/ij^gN, such that 

Xj £ a{M nj ), fij e a{L nj ), 

and Xj — > fi, (ij — > fi as j ' oo. 

A consequence of this theorem is that if a subsequence of eigenvalues of the regularized operators 
M n converges to some point \x which is not an eigenvalue (spectral inexactness) then the L n will 
also possess a subsequence of eigenvalues converging to the same point fi. Moreover, if subsequences 
of eigenvalues of M n and of L n converge to the same point then at least one of the subsequences 
is spectrally inexact, because the boundary conditions ( |l4| ) and (|l^) ensure that M and L have no 



shared eigenvalues. Theorem 3.1 therefore gives us a test for spectral exactness: if only the M n , and 



not the L n , possess eigenvalues accumulating at /x, then /i must be an eigenvalue of M. 



Proof of Theorem 3.1. Let A be any sufficiently small annulus surrounding /i and let T be a closed 
contour in A surrounding fi. For any function / which is meromorphic in a simply connected open 
set containing T we denote by Nz(f,T) the number of zeros of / inside T and by Np(f,T) the 
number of poles of / inside T. Rouche's Theorem gives 

N z (m,T) - N P (m,r) = ^ f ^^dA. 
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In view of the identity m(X)£(X) = —1 we have Nz(m,T) = Np(£,T) and 1/m = —£, so we can 
write this 

N P (i, T) - N P {m, T) = — — [ m'(X)£(X)dX. (19) 

2m J r 

As /i does not lie in the spectrum of M, we have Np(m,T) = for all sufficiently small annuli A. 
Hence 

Np(e,T) = -— [ m'(\)l(X)d\. (20) 



2m _ , 

By arguments similar to those which gave ( fL9| ) we have 

iVp(C,r) -iVp(m n) r) = J^m! n (\)£ n (X)dX. 

The uniform convergence m„ — > m implies uniform convergence of m' n to m' (by the Cauchy integral 
representation of the derivative) . Combined with the uniform convergence £ n — > £ this yields 



N P (£ n ,T) -N P {m n ,T) = --L / m'(\)£(\)d\ 



for all sufficiently large n. Combining this with ( pOj ) we have, for all sufficiently large n, 

N P (i,T) = N P (£ n ,T) - N P (m n ,T). (21) 

In the case of possibility (a) , the sequence M n does not have any eigenvalues converging spuriously to 
which is a non-eigenvalue of M. Since eigenvalues of M n are poles of m n we have Np(m ni T) = 
and (]2l]) then shows that the sequence L n is spectrally exact for L near /i. Thus we have spectral 
exactness near /i for both M and L. In the case that (a) is not true, then for some arbitrarily small 
annuli A and arbitrarily large n £ N we will have 

Np(m n , T) = no. of eigenvalues of M n inside T > 0, 

and, from @, 

N P (£ n ,T) -N P (m n ,T) = N P (£,T) > 0, 

which shows that Np(£ n ,T), the number of eigenvalues of L n inside T, is at least 1. Thus L n also 
has an eigenvalue close to /i. This gives possibility (b). □ 

3.2 Conditions for m-function convergence 



We now examine the hypotheses ra n — > m and £ n — > £ of Theorem 3.1. Under what conditions do 
these hold? 

We consider first Sims Cases II and III. The following result explains how to choose the boundary 
condition (|l^) to ensure that ([11]) and ([l8|) hold. 

Lemma 3.2 Suppose that the differential equation is of Sims Case II or III. Let X' G A ?) jf be 
fixed. Express the boundary conditions at x = b for M in terms of an L^-solution ip(-,X') = 
6{-, A') + m(A')</>(-, A') of the differential equation A4ip — X'?p, in the form 

[y,i;(;X , )}(b)=0, 

^ eq. (4- 11)] where [•, •] denotes the usual Wronskian [u,v] := p(uv' — u'v). Then appropriate 
boundary conditions |m) are given by choosing 



(cos/3„,sin/3„) = const. (pip' (b n , A'), -ip(b n , X')) (22) 
so that (|7qj is simply the condition [y, ip](b n ) = 0. 
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Proof Define tp n (-, A') = 6{-, A') + m n (X')(f>(- , A'), in which m n is chosen so that t/; n satisfies ( plf ) 
with /3„ given by (p^). From the definition of the f3 n , the fact that ip n satisfies ( |l6| ) may be written 

as 

[Va-,A'),V(-,A')](M = 0. 
Substituting Vn('> A') = 9(-, A') + m n (X')(f)(- , A') into this equation yields 

[fl(-,A'),V»(-,A')](6 n ) 

m " (A) = l^n^A0](M- (23) 

In the identity [^(-, X'),ip(; X')](b) = 0, replace the first instance of A') by 9(-, X') + m(X')4>(-, A'), 
and hence obtain 

m(A) = -[0(-,n^(-,AO](6)- (24) 

Comparing (§|) with (§|) establishes (0) when A = A': 

lim m„(A') = m(A'). (25) 

n — >oo 

Now Brown et al. j|, Corollary 3.4, eqn. (3.4)] give a formula which allows us to extend this result 
to other values of A: 

m(A') - (A - A') f b w(x)9(x, XU(x, X')dx 

m(A) = — - — [ Ja — V — ' . (26) 

1 + (A - A') f° w(x)(l>(x, X)ip(x, X')dx 

This formula possesses the regular-interval analogue 

m n {X') - (A - AQ / a b " W (x)g(x, A)V>„ (x, X')dx 

m n (X) = 7 . (27) 

1 + (A - A') J a " w(x)<f>(x, X)ip„ (x, X')dx 

These formulae hold at any point which is not an eigenvalue of M or of M n , respectively. Moreover, 
since the equation is in Sims Case II or III, all of 6(-, A), 9(-, A'), </>(■, A) and </>(-, A') lie in L^,[a, b). 
Using ip n (-,X') - i>(;X') = (m„(A') - m(X'))<p(; A') it follows from @ that ip n (-,X') -> ^(-,A') in 
L^Jfl, b). Hence, combining ( p6| ) and (|27|), we obtain the convergence 

lim m n {X) = m(X) 



at any point A which is not an eigenvalue of M. This establishes (17), and (|18|) is proved similarly. 
□ 



Of course, Theorem 3.1 requires more than just pointwise convergence, and so it is fortunate 
that the following stronger result holds. 

Theorem 3.3 Suppose that the problem is Sims Case II or Sims Case III at x = b and let the 
hypotheses of Lemma 3.i hold. Then for A in any compact set K, C C not containing eigenvalues of 
M, 

lim m„(A) = m(A), (28) 

the convergence being uniform over IC. 

Proof That m n (A) — > m(A) pointwise on JC has already been proved in Lemma [T^. The uniformity 
of the convergence depends on having a uniform bound on the L^, norms of 6(-, A) and (/)(■, X) for 
A G K. This can be obtained by a standard variation of parameters argument, expressing the 
solutions in terms of 9(-,X') and (p(-,X'): see Sims |l6|, section 3, Theorem 2] and also [|], Remark 
2.2]. □ 
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For the functions £ n and £, a result exactly analogous to Theorem [}7| is clearly valid, the only 
difference being now that fC must not contain eigenvalues of L. 

We turn now to Sims Case I. In order to handle this case it is necessary to know more about the 
behaviour of the solutions of the differential equation. Suppose that for some A € C, the differential 
equation possesses 'small' and 'large' solutions. We shall assume that the small solution is the 
(unique up to scalar multiples) square integrable solution tp(x, A), and we denote the non-unique 
large solution by T(x, A). By 'small' and 'large' we mean that these solutions satisfy the condition 

lim^4=0- (29) 

x^b I [X, X) 

Clearly T is not unique: T + ip, for example, is also a 'large' solution in the sense of (p9|). The 
solutions 9 and (f> of (|7|) can clearly be written in terms of tp and T: 



9(x, A) = c\ijj{x, X) + c 2 T(x, A), 
</>(x, A) = di^(x, A) + d 2 T (x, A), 



(30) 



in which the constants c\ , C2 , c?i and d 2 are given by 

ci = ((cosa)pT'(a, A) - (sina)T(a, A)) /W, c 2 = (-(cosa)p^'(a, A) + (sin a) ip (a, A)) /W, (31) 
d x = ((sina)pT'(a, A) + (cos a) T (a, A)) /W, d 2 = (-(sin a)ptp'(a, X) - (cos a)ip(a,X)) /W, (32) 

where W — p( - 0T / — ^>'T) is the usual Wronskian. Suppose that m n is defined by the requirment 
that the solution 

V>„(-,A)=0(-,A)+m n (A)<K-,A) 
satisfy the boundary condition ip n (b n , A) = 0. Then 

m «(A) = -— - — --. (33) 

<P\P n , A) 

Now combining (|29| ) with ( ^p| ) we have 

0(b n , A) ~ c 2 T(6„, A), <j>(b n , A) ~ d 2 T(6„, A) 
for large n. Combining this with (^) yields 

m n {\) ~ 

d 2 

for large n. Together with ( |3l| ) and (|3^) this yields, for large n, 

-(cos ai)pi>' (a, X) + (sma)ip(a,X) 

Tn n (A) ~ — : : ; rr ; TT7 7^ = m (X), (34) 

{ ' (sin a)pip'(a, A) + (cos a)ip(a, A) V y ' V ; 

the last equality in ([mJ) being an immediate consequence of Definition 4.10]. From these consid- 
erations the following result is clearly true. 

Theorem 3.4 Suppose that the differential equation is of Sims Case I type at x — b. Let K, C C be 
a compact set such that for A € £ the square-integrable solution ip( x i A) of the differential equation 
exists and is an analytic function of X. Suppose moreover that there exists a second solution T(x, X) 
such that 

ip(x,X) 

x-^b 1 (X, X) 

the limit being uniform with respect to X G /C. Suppose that the domains of the operators M n are 
determined by the Dirichlet conditions y(b n ) = 0. Then we have the convergence 

lim m n (X) = m(X) 

n — too 

uniformly for X € tC. 
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Clearly a similar result holds for the functions £„(X) and their convergence to the function £(X). 



Remark The result of Theorem 3.4 will also hold if the domains of the M n are defined by certain 



other boundary conditions at x = b n . Suppose that a boundary condition 

y(b n ) cos/3„ + y'(b n ) sin/3„ = 

is imposed, where the (3 n are complex numbers. Then it may be shown that the result continues to 
hold provided 

ip(b n , A)cos/3„ +-0'(&n,A)sin/3„ 
™^oo T(6„,A) cos/3„ + T'(6„, A)sm/3„ 

locally uniformly with respect to A. In problems where ip'(x, A)/T'(x, A) — > as x — > b, one would 
have to choose the [3 n quite carefully for (^) to fail. 

3.3 A simple test for spectral inclusion 

In the selfadjoint case, spectral inclusion is usually very easy to prove. In fact, suppose T is a 
sclfadjoint operator on a domain D(T) in a Hilbert space H and let (T n ) be a sequence of operators 
with domains (D(T n )) which converge pointwise to T on some set C C n m£ N U n > m D(T n ): 

lim \\T n f - Tf\\ =0 V/ G C. (37) 

n — >oo 

Then provided C is a core of T - in other words, provided the set (C,TC) is dense in the graph of 
T - the sequence (T„) will be spectrally inclusive for T: every eigenvalue of T will be the limit of 
some sequence (A 1 -™') in which X^ n > lies in the spectrum of T n . 

In the non-selfadjoint case a result of such generality does not seem to exist, although some of 
the results of Harrabi come quite close. We shall examine some corollaries of Harrabi's work, as 
well as a standard result from Kato Jl3|], in section ^ below. In this subsection, however, we shall 
show that spectral inclusion always h olds in Sims Cases II and III, and in Sims Case I in those parts 
of the complex plane where Theorem p.4| holds. 



Theorem 3.5 (Test for Spectral Inclusion) In the notation of Theorem 3.1, suppose that /igC is 
an isolated eigenvalue of M. Suppose also that m n (A) — ► m(A) as n — * oo uniformly on any compact 
annulus surrounding \x. Then there exists a sequence (A*-"') in which A*-™' lies in the spectrum of 
M n , such that lim^oo A^ 1 ' = fi. 

Proof Given e > sufficiently small, surround /j, by an annulus A whose outer radius is at most e 
and let T be a circular contour surrounding fi and contained in A. Since /i is a pole of to, we may 
assume by taking e sufficiently small that A contains no zeros of to. The uniform convergence of 
m n to to on A then guarantees that for all sufficiently large n, m n is bounded away from zero in A, 
so l/m n also converges uniformly to 1/to in A. Moreover, Cauchy's integral representation of the 
derivative implies that m n converges uniformly to ml on A. Hence we have 

1 f <(A) „ 1 f ro'(A) , , 



2ni J T TO n (A) 2iri J r to(A) 

and the fact that both sides of this equation are integers gives 



2iri J r to„(A) 27ri J T to(A) 

for all sufficiently large n. 



1 f<iX) dX = ±[^ d X (38) 
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From the argument principle, if we let v denote the algebraic multiplicity of /i as an eigenvalue 
- and hence as pole of m - we have 



1 f m '^-dX. (39) 



2iri J r m(A) 

Combining ( pH ) and (^) shows that M n also has eigenvalues of total algebraic multiplicity v inside 
the contour F, for all sufficiently large n. □ 
In section |5| we shall give some examples in which the Eastham-Levinson asymptotics |8| allow 
us to verify the hypothesis (^5|) and hence apply the test for spectral exactness given in Theorems 



3.1 and 



4 Pseudospectral inclusion and spectrum-of-sequence inclu- 
sion 

In this section we consider a sequence (T ra ) of operators on a Hilbert space H . Let T be some other 
operator on H. We denote by D(T n ) the domain of T n and by D(T) the domain of T. We shall be 
interested in two different types of convergence of T n to T: strong convergence, in which T n f — > Tf 
for each fixed /, and norm resolvent convergence, in which || (XI — T„) _1 — (XI — T) -1 1 — > 0. Strong 
convergence is usually observed when a problem of Sims Case I is regularized by a sequence of 
interval truncations, while the stronger property of norm resolvent convergence is observed when 
the problem being regularized is of Sims Case II or III. Strong convergence generally results in a very 
weak type of spectral approximation which is given by Theorem [L3 below; in practice, for differential 



equation eigenvalue problems, this is probably not as useful a result as Theorem 3.5. Norm resolvent 
convergence, on the other hand, gives a spectral exactness result (Theorem 4.5) which Theorems 
3.1 and 3.3 do not give: Theorems 3.1 and 3.3 do not preclude spectral inexactness in Sims Cases 
II and III, they merely give a test for spectral inexactness, whereas Theorem 4.5 precludes spectral 
inexactness. 

The following definition is standard. 

Definition 4.1 A set C C D(T) is called a core ofT if, for every x G D(T) and e > 0, there exists 
i f GC such that 

\\x — x e \\ < e, \\Tx — Tx e \\ < e. 
The following definition is also required. 

Definition 4.2 The spectrum of the sequence (T n ), denoted a({T n }), is the set 

<i({T n }) = {X e C \ lim WiXI-T^W =+co}. 

n — >oo 

Note that for selfadjoint operators, cr({T„}) can contain only points which are limit points of se- 
quences of the form (A 1 -™-*) in which X^ lies in the spectrum a(T n ) of T n : 



a ({Tn}) C n me NU„> m cr(T„). 



(40) 



To see this, let A € a({T n }) and let a n := \\(XI — T n ) _1 |j _1 , so that a n — > as n — > oo. By the 
spectral calculus for the selfadjoint operator T„, given e > there certainly exists a point of a(T n ) 
within distance a n + e of A: in particular, choosing e = — we can find A< n ) e o-(T n ) such that 
| A — A n | < a n + —. Hence for any integer m, 



A G U n > m a(T n ). 
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This proves ©. For the non-selfadjoint case the spectral calculus no longer holds (unless the T n 
happen to be normal) . The following result - a simple modification of a result of Harrabi |y] - thus 
provides for non-selfadjoint operators as close an analogue of the spectral inclusion result of Reed 
and Simon jlfj, Theorem VIII. 24] as is possible, in general. 

Theorem 4.3 Let C be a core ofT and suppose that 

C c n mG N U„> m D(T n ), 

so that if f G C then T n f is defined for all sufficiently large n. Let 

o~ a {T) := {A £ C | 3(a;„)„ g N C H such that \\x n \\ = 1 Vn and lim ||(AI — T)x n || = 0}. 

n — >oo 

Suppose that for each f G C we have linin^oo \\T n f — Tf\\ = 0. Then 

<Ja(T)Ca({T n }). (41) 

Proof Suppose that A does not lie in a({T n }). Then there exists M G R + and a monotone increasing 
sequence (rij)j^ such that 

UXI-Tn^WKM. 
Now let feC. Then f = (XI - T nj )^ 1 (XI - T nj )f, whence 

ll/H < M|| (XI - T nj )f\\ -» M||(A7 - T)f\\ as j -> oo. 

Since this holds for all / G C, it follows from Definition |4.1| that 

ll/ll <M||(A/-T)/|| VfED(T). 

Hence A does not lie in a a (T), which proves the result. □ 
Because the result ( ^0| ) does not hold, in general, for non-selfadjoint operators, it is not generally 
possible to replace a({T n }) in ( fil| ) by r\ me fqL) n > m a(T n ) . However, the following result concerning 
pseudospectra can be proved. 

Theorem 4.4 Let e > 5 > 0. Let 

a e (T„):={AeC| || (XI - T n )- X || > e" 1 }, (42) 

as(T) := {A G C | ||(A/ - T) _1 || > cT 1 }. (43) 
Let C be a core ofT satisfying the same hypotheses as in Theorem and suppose that for all f G C 



we have linin^oo \\T n f — T/|| = 0. Then 

a s (T) Climmfo-e(T n ), (44) 

n — >oo 

where 

liminf a e (T n ) := U mG N 

n — »oo 

Proof We start by defining 

a s ({T n }) := {A G C | liminf IKAI-T,,)- 1 1| ^r 1 }. (45) 

n — >oo 

Suppose that A does not lie in crs({T n }). Then there exists 7 G (0, S^ 1 ) and a monotone increasing 
sequence (rij)jeN such that 

ll(A/-T„ j )- 1 ||< 7 <r 1 . 
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Let feC, and write f = (XI - T^y^XI - T n .)/, giving 

11/11 <7ll(A/-T n .)/||. 

Letting j — > oo gives 

11/11 <7ll(A/-T)/|| (46) 

Eqn. holds for all f e C and hence, since C is a core of T, for all / e D(T). In particular this 
implies that XI — T is invertible and 

\\(XI-T)- l f\\< 7 \\f\\ 

for all / e D(T). This implies that ||(A7 - T) _1 | < 7 < 6' 1 , and so A does not lie in u s {T). We 
have thus proved 

a s (T)Ca 5 {{T n }). (47) 
The result will be proved if we can show that for S < e, 

<Ts({T n }) C liminf a e (T n ). (48) 

n — *-oo 

To do this, suppose that fi does not lie in liminf n _>oo a e (T n ). By definition, 

liminf a e (T„) = U meN n„> m {A G C | ||(AJ - T n ) _1 | > e" 1 }, 

n— »oo 

and so for each m G N, /x does not lie in 

n n > m {AeC| IKAi-T,,)- 1 !! ^e" 1 }. 

In other words, there exists a subsequence (T nj ) jgN such that 

IKAtT-T^)- 1 !! <e- 1 J jGN. 

Hence by definition, 

liminf || ((/X/-T,,)- 1 1| < e" 1 < *~ x . 

n — >oo 

From ( ^5| ) we have clearly proved /it does not lie in as({T n }). This establishes (ff8|), and our proof is 
complete. □ 

Theorem 4.5 Let z e C &e /raed. Suppose that \\(zl - T n )- X - (zi - T) _1 || -» as n — > 00. Then 

limsupcr(T„) C ct(T), (49) 

n — >-oo 

where limsup is defined by 



limsupa(T n ) = n, meN U n > m a(T n ). 
n— >oo 

Proof Let i?„ = (zJ - T„)" 1 and let R= (zl- T)~ x . Then ||i? n - -> as n -> 00. From the 
results in Kato g IV, §3, p. 208] it follows that 

limsup o-{Rn) C o-(ij). (50) 

n — >oo 

However the spectrum of i? is related to the spectrum of T by 

A e a{T) if and only if (z - A)" 1 G a{R). 
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[N.B. for our applications, oo will be an accummulation point of cr(T) and so will lie in a(R).] A 
similar relationship holds between a(T n ) and cr(i?„). Thus (pQ) implies □ 



We shall now show that in Sims Cases II and III, the hypotheses of Theorem 4.5 are satisfied by 
taking T n = M n and T = M, where M n and M are the operators of Theorem 3.1. Thus Theorem 



4.5 will supercede Theorem 3.1 in Sims Cases II and III as a guarantee that spectral inexactness is 



impossible provided the boundary conditions are correct. Spectral inclusion still holds by Theorem 
|3.5[ Combining all these results will show that in Sims Cases II and III, provided we generate the 
M n using the boundary conditions described for these cases in Lemma 3^ , we have spectral inclusion 
and spectral exactness (Theorem [4.7| below). 



Theorem 4.6 Consider a differential expression of Sims Case II or Case III type at x = b. Let M 
be a realization of this expression through analytic continuation of an m-function which is a limit of 
functions m n (-) for realizations M n of the differential operator defined on intervals [a, b n ], b n f b as 
n y oo, as described in Lemma 3A. Let ip n (x,X) be the solution of the differential equation defined 
by 

ip n (x, A) = 6(x, A) + m n (X)(f>(x, A), 
and let the G n (x, y, A) be the Green's functions given by 



G„(x,y,X) = 



-4>{x,X)ip n (y,\), a<x<y<b, 
-ip n (x,X)<j)(y,X), a<y<x<b. 



Let Rn(X) be the extension to L^[a, b) of (XI — M n ) 1 defined by 



(Rn(X)f)(x) 



G n (x,y,X)f(y)w(y)dy, f £ L w [a,b). 



(see eqn. (4-2)]). Fix z £ C and suppose that z does not lie in the spectrum of M or of any of 
the M n for sufficiently large n. Let R = (zl — M) . Then \\R(z) — R n (z)\\ — > as n — > oo. 



Proof From j|, eqn. (4.2)] we know that 

(R(z)f)(x) -- 



G{x,y 1 z)f(y)w(y)dy, 



in which ip(x, z) 



[x, z) + m(z)<f>(x, z). Thus 

r-b 



G{x,y,z)f(y)w(y)dy+ / (G(x,y,z)-G n (x,y,z))f(y)w(y)dy. 



(R(z)-R n (z))f)(x) 



Because the differential equation is of Sims Case II or Case III, we know that both #(-, z) and </>(•, z) 
lie in L^[a, 6). This implies that 



w(x)dx I w(y)G(x,y, z) dy < +oo, w(x)dx / w(y)G n (x,y, z) dy < +oo 



In particular, therefore, 



lim I w(x)dx I \G(x,y,z)\ w(y)dy = 0. 



The bound 



b„ 



\\R{z) ~ R n {z)\\ 2 < 2 I w{x)dx JJG{x,y,z)\ 2 w{y)dy 



2 w(x)dx w(y)dy\G(x,y,z) - G n (x,y,z)\' 
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now yields 



lim \\R(z) - Rn(z)\\ <2 lim w(x)dx w(y)dy\G(x,y,z) - G n (x,y,z)\" 

n — >oo n — >oo ' 



-(m(z) - m n (z))4>(x, z)(j){y, z), a < x < y < b, 
-(m(z) - m n (z))9(y, z)<j)(y, z), a < y < x < b. 



We now use the formulae 

G(x,y,z) - G n (x,y,z) = 
to obtain 

lim || J R(z)-i?„(z)|| 2 <2(lim \m n (z) - m(z)\ 2 ) I fw(OHL z)\ 2 dA I [ w(0\cf>(Z, z)\ 2 d£ 
Since linin^oo \m n {z) — m(z)\ = 0, this proves the result. 



□ 



Theorem 4.7 (Spectral Inclusion and Spectral Exactness for Sims Cases II and III). Let M n and 
M be as in Theorem \j.b\ Then 

(a) for every A in the spectrum of M , there exists a convergent sequence (X^) ne fi, with \^ in the 

spectrum of M n , whose limit is A; 

(b) if (A' n )) n eH is a convergent sequence with limit A and A*™-* lies in the spectrum of M n for each 
n, then A lies in the spectrum of M. 



Proof By Theorem 4.6, the hypotheses of Theorem 4.5 a re sa tisfied. This immediately gives (b). 
Tur ning to (a), we observe tha t th e hypotheses of Theorem 3.3 are satisfied. The result of Theorem 
3.3 allows us to use Theorem 3.4, which in turn allows us to use Theorem 3.E. The conclusion of 
Theorem 3.5 is precisely (a). □ 



5 Examples 

We illustrate the results of the preceding sections with some numerical examples. 

Example 1 An equation of the form 

-y" + c 2 y = \w(x)y, x € [0, oo), 

in which 5R(c) ^ and w(x) = exp(— 3|5ft(c) \x), is easily checked to be in Sims Case II at infinity. 
Letting v(x) = exp(— |5R(c)|x) we can define an operator M by (My)(x) = w(x)~ 1 {— y" + c 2 y} 
for y £ D(M), where the boundary conditions for D{M) are 

2/(0) = 0, [y,w](oo) = 0. 

The corresponding operator L has domain D(L) specified by the boundary conditions 

?/(0) =0, [y,»](oo)=0. 

For the operators M n and L n on finite intervals [0, b n ] the boundary conditions at the origin 
will be the same as for M and L respectively, while the boundary condition (^2|) at x = b n 
will be given, according to Lemma 3.2, by 

[y,v](bn) = 0. 
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Using the code described in [|10[ we computed the eigenvalues of the operators M n and L n in 
the box with corners 100, 100(1 + i), lOOi, 0. The results, shown in Table 1, indicate that 
the eigenvalues of the M n and of the L n converge to distinct points in the box, and so our 
test for spectral exactness suggests that the operator M has eigenvalues close to 94.4890 + 
£28.8595 and 24.21335 + £14.11108, while L has eigenvalues close to 3.1163595 + £5. 808222 and 
51.51888 + £21.3277. Of course, this is what we would expect for a Sims Case II problem, by 



Theorem 4.7 



Note that these eigenvalue problems can be formulated as compact perturbations of selfadjoint 
eigenvalue problems, although it is not immediately clear how one might use this to obtain 
spectral inclusion and/or exactness results. 
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K 


Eigenvalues of M n 


Eigenvalues of L n 


1 


5 


24.21311+il4.10915 
94.4880+i28.8342 


3.1163619+i5.808222 
51.51912+i21.3277 


2 


10 


24.21334+il4.11103 
94.4891+i28.8584 


3.1163595+i5.808219 
51.51879+121.3274 


3 


15 


24.21333+il4.11106 
94.4888+i28.8590 


3.1163591+i5.808221 
51.51875+121.3276 


4 


20 


24.21335+il4.11108 
94.4890+i28.8595 


3.1163595+i5.808222 
51.51888+i21.3277 



Table 1: Example 1 on intervals [0,6„], using code of |[o) with TOL = 10 



Example 2 We consider the (now rather infamous) rotated harmonic oscillator problem 

-y" + c 2 x 2 y = Xy, a;e[0,oo), 3/(0) = 0, c e C, ffi(c) > 0, 

(see Davies (6|). This problem is of Sims Case I at infinity and its eigenvalues are given by 

A fe = c(4fc + 3) fc = 0,l,2, 

It is known that the higher index eigenvalues are very ill-conditioned (when c 2 is not posi- 
tive). Denoting by M the operator associated with this problem, this ill-conditioning may be 
explained by the fact that \\(M — z/) _1 || is extremely large in very large neighbourhoods of 
these eigenvalues, making it numerically difficult to determine the precise location of the poles 
of || (M — z/) _1 ||, which are the eigenvalues. 

are satisfied for this 



On the other hand, it is easy to verify that the hypotheses of Theorem 3.4 



problem, so Theorem [ST] still gives a valid test for spectral inexactness. 

One might expect that this would be of rather academic interest, given that the operators 
M n and L n on the truncated intervals [0, b n ] will themselves have very ill-conditioned higher 
index eigenvalues. To some extent this is correct. However, in Table [| we show the result of 
truncating the interval to [0, 20] and locating all the eigenvalues in a rectangle in the complex 
plane with bottom right-hand corner A = 100 and top left-hand corner A = 90z. The boundary 
conditions used were y(0) = = 3/(20) for the M n problem and y'{0) = = 3/(20) for the 
L n problem. The spurious eigenvalues are marked with asterisks. One can see quite clearly 
that these eigenvalues distinguish themselves by being almost invariant under the change of 
boundary condition at the origin. Indeed, in all but one case the relative differences are less 
than the tolerance which was used in the computations (10~ 5 ). There is also a problem with 
'missing' eigenvalues in this table: M n ought to have an eigenvalue close to 73 + £50 and 
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L n ought to have an eigenvalue close to 77 4- i55, both of which are missing. Thus the ill- 
conditioning of these problems may induce spectral inexactness, for which we seem to be able 
to test by changing the boundary conditions, but it can also cause a lack of spectral inclusion, 
which is rather more difficult to spot. 



Eigenvalues of M n 



Eigenvalues of L n 



4.3278454 

10.098296 

15.868726 

21.639117 

27.409413 

33.179628 

38.949823 

44.720108 

50.490469 

56.260702 

62.032125 

67.763082 

70.791811 - 

76.994286 

72.268485 - 

73.809759 - 

75.474904 - 

87.360734 

98.348465 



i 3.1193175 

- i 7.2784044 

- i 11.437476 

- i 15.596455 

- i 19.755403 

- i 23.914366 

- i 28.073359 
i 32.232282 
i 36.391057 
i 40.549373 

- i 44.715888 

- i 48.620308 
i 54.525167* 

- i 49.224379 
i 64.384873* 
i 74.921450* 
i 86.054406* 

- i 47.089232 

- i 44.849646 



1.4426265 
7.2130633 
12.983488 
18.753885 
24.524267 
30.294667 
36.065062 
41.835242 
47.605269 
53.375120 
59.146275 
64.868061 
70.836011 - 
71.821431 
72.268524 - 
73.809759 - 
75.474905 - 
82.109895 
92.771537 



- i 1.0397661 

- i 5.1988778 

- i 9.3579536 

- i 13.517014 

- i 17.676070 

- i 21.835034 

- i 25.993830 

- i 30.152542 

- i 34.311200 

- i 38.469967 

- i 42.628491 

- i 46.792270 
i 54.242901* 

- i 49.862180 
i 64.384384* 
i 74.921449* 
i 86.054405* 

- i 48.176485 

- i 45.977747 



Table 2: Testing for spectral inexactness on the rotated harmonic oscillator 



Example 3 Consider the problem of locating resonances of the equation 

- y" + 16x 2 exp(-x)y = Ay, y'(Q) = 0, ate[0,oo). 

Using the 'complex scaling' method, the resonances of this problem are of the form e" 
where \xq is an eigenvalue of the non-self adjoint problem 

- z" + 16x 2 e 2 ' 9 exp(-xe i9 )z = fj,z, z'(0)=0, x€[0,oc), 



(51) 



Me 



(52) 



(see Hislop and Segal [f[2|). The rotation angle 6 > must be such that the function x i— > 
16x 2 e 2ie exp(— xe l8 ) lies in L 1 [0, oo), and in particular therefore 9 < n/2. Resonances have 
the property that e~ 2te fj,g is independent of 8, so in general not all eigenvalues of (|5^) yield 
resonances: one should carry out the computations for at least two different values of 6 to 
identify resonances. 



In addition to the complications caused by the fact that some eigenvalues of (52) do not 
correspond to resonances, we have the additional problem that (52) is a singular problem and 
must be regularized by interval truncation. This truncation process might introduce spurious 
eigenvalu es, not corresponding to eigenvalues of (|52|) , which we need to be able to detect. 
Theorem 3T gives a way to do this. 

with boundary condition 
with y'(0) = and the 



Using a rotation angle — 1.1 and a truncated interval [0,100 
?/(100) = we computed both the eigenvalues of the equation in 
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eigenvalues for the same equation but with the boundary condition y(0) = using the code 
of pJ]. We asked the code to find, in the /i plane, all the eigenvalues in the box with corners 
(-0.01,0.01), (-0.01,5), (-10,5), (-10,0.01), with a tolerance of 10~ 6 . The results, rotated 
back into the A plane via A = e~ 2t9 (i, are shown in Table 0. 



Alleged resonances 
with y(0) = = ?/(100) 
on [0, 100] 


Alleged resonances 
with y'(0) = = 2/(100) 
on [0, 100] 


2.429823932+i 2.95502902 
3.869964809-i 0.74439879 
0.554661821+i 0.66915540 


2.429823937 +i 2.95502903 
3.869964804 -i 0.74439879 

0.554661961 +i 0.66915556 
2.861786706 -i 1.6 x 10" 6 



Table 3: Testing for spurious resonances due to interval truncation 



Of the four alleged resonances found with y'(0) = 0, three are virtually unchanged when the 
boundary condition is changed to y(0) = 0. Theorem 3.1 indicates that these are probably 
spurious. This is obvious for the two which have positive imaginary parts, as resonances lie 
in the lower half plane by definition; however, without Theorem 3.1 it would not have been 
obvious for the alleged resonance at 3.8699648 — iO. 7443988. For this problem we believe that 
the only genuine resonance found, for boundary condition y'(0) = 0, is the one at 2.861786706 — 
il.6 x 10 -6 . In fact, of the four alleged resonances this is the only one which is invariant under 
a change of the rotation angle 9; however, in general it is not clear that a spurious resonance 
generated by interval truncation would always fail to be invariant under change of 8. 
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